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Detection of topological transitions by transport through molecules and nanodevices 
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We analyze the phase transitions of an interacting electronic system weakly coupled to free- 
electron leads by considering its zero-bias conductance. This is expressed in terms of two effective 
impurity models for the cases with and without spin degeneracy. We demonstrate using the half- 
filled ionic Hubbard ring that the weight of the first conductance peak as a function of external flux 
or of the difference in gate voltages between even and odd sites allows one to identify the topological 
charge transition between a correlated insulator and a band insulator. 
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Progress in nanotechnology has made it possible to 
perform transport experiments on systems as small as 
single molecules |lj. Metallic (2| or semiconducting 
[2|,y,y| quantum dots (QDs) can now be assembled into 
artificial molecules |4|| or solids l2| . QD molecules of dif- 
ferent materials and sizes are now being investigated and 
a wide range of new QD systems is expected to be synthe- 
sized in the near future ?, 4]. The transport properties 
of a finite chain of 15 QDs were first measured over 10 
years ago S] , and the metal- insulator transition has been 
studied experimentally in a hexagonal lattice of Ag QDs 
[21 • These advances open the route for new approaches 
to investigate novel phenomena and theoretical concepts 
in interacting electron systems. 

Here we focus on one such phenomenon, the topological 
phase transition associated with a parity change of the 
ground state. The ionic Hubbard model (IHM) with al- 
ternating diagonal energy ±0^ ^^^ received much recent 
attention 0, S B B ll^ ■ It was proposed to describe the 
neutral-ionic transition in organic charge-transfer salts 
[llj, and later applied to model ferroelectricity in per- 
ovskites [g], [ij]- At half filling and in the atomic limit 
(t — > 0), the ground state is a band insulator (BI) for 
U < A, but is a Mott insulator (MI) for U > A. In 
one dimension, with non-zero hopping t, a spontaneously 
dimerized insulator (SDI) phase appears between BI and 
MI. With increasing U, one finds first a charge transi- 
tion sX U = Uc from BI to SDI, followed by the closing 
of the spin gap at Us > Uc, where the transition to the 
MI occurs. Although in finite systems conventional or- 
der parameters such as charge and spin structure factors 
vary continuously at a phase transition, in a system of 
L sites one may define charge and spin topological num- 
bers which change discontinuously at Uc{L) and Us{L) 
|g|. The charge Berry phase has a step at Uc{L), where a 
parity change of the ground state occurs for periodic (an- 
tiperiodic) boundary conditions if L ^ Am {L = Am + 2, 
m integer) [a,l3- 

Here we show that this charge transition can be de- 




FIG. 1: Interacting electron system on a flux-threaded ring, 
connected by weak links (f') to two conducting leads. 



tected using the total intensity or width w of the first 
peak in a zero-bias conductance measurement performed 
on a flux-threaded ring. The importance of the ring ge- 
ometry is that the topological transition is absent in an 
open system |3- In one of the two parity sectors, which is 
selected using an applied gate voltage, w ^ for applied 
flux (j) —>■ if L = Am [(j) — > (j)o/2, where ^o = hc/e is 
the flux quantum, if L = 4to -|- 2). Thus the transition, 
which is observed by varying parameters such as A, may 
also be studied in ring-shaped molecules where it is not 
possible to attain a significant threading fiux. 

The configuration for the proposed experiment is il- 
lustrated in Fig. 1. The interacting electron system on 
the ring is connected to conducting leads through two 
sites, labeled and n, by weak links with hopping i', 
and either two materials or two different gate voltages 
at alternating sites may be used to model the IHM. It 
is essential to distinguish between two cases depending 
on the spin degeneracy of the isolated interacting system 
for odd particle number. We thus perform two different 
calculations of the conductance. The first assumes that 
spin degeneracy is lifted by a Zeeman interaction, but 
allows orbital degeneracy of the states, which is impor- 
tant when including interference effects arising in a ring 
geometry [13. In the second, where spin degeneracy is 
retained, we map the problem into an effective Anderson 
model (EAM) and express the conductance as a function 
of the spectral density of the latter. 

If the ground state \g) of the interacting system is non- 



degenerate, for small t' the links can be eliminated by a 
canonical transformation, leading to an effective, non- 
interacting Hamiltonian for the two leads in which the 
two "impurity" sites i connected to the links (—1 and 1 
in Fig. 1) have an energy shift Aei(w), and are connected 
by an effective hopping ieff(w), 

HcS = ^ efc^c{,^^Cfe^cr + Ae_ic^ig.c_io- + Aeicl^cicr 



kR.a 



efc 



fl4fl<TCfeH<T + ^{toSc\„C^la + H.C.). (1) 



L (R) refers to the lead containing the site —1 (+1), and 
the impurity parameters may be expressed in terms of the 
Green functions for the isolated ring gij{uj) — {{cia', Cj^)), 



Ae_i(t. 



ieff (w) = t'^9no{uj). (2) 



The transmittance of Hcs may be calculated in dif- 
ferent ways. Without affecting the essential results, we 
assume identical leads and, for their states without sites 
(—1,1), we consider two models for the density of states 
p{lu) and the hybridization Vj{lu) with sites (—1,1). I) 
p{llj) and Vj{llj) constant: for j ~ ±1 and teff — 0, 
g]f{uj) = l/iuj-Ae-i + iT) with T = irpV^. II) Leads de- 
scribed by a tight-binding model with nearest-neighbor 
hopping t, where 5|f°(w) = 1/(^/2 - Ae^ + iT{uj)) with 
r(a;) — 7rp{LLi)V-^{uj) — y^f^ — ^"^ j \. By introducing an 
integer ttt, — 1 or 2 for cases I and II, 



T^'^.yu) 



(3) 



4r2(a;)|i,ff| 



f^-A. 
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zr(u;))(^-Aei+ir(a;))-|ieffPp- 



K) is a gate voltage which changes the on-site energy of 
all sites of the ring by —eVg. In case II this equation is 
exact (Vi') in the non-interacting system, and generalizes 
a previous result [l3- Eq. © also generalizes previous 
approaches in which only one intermediate state of the 
ring is included [ij, Il5l. a nd is appropriate for the study 
of interference effects |l3l | . Spin degeneracy limits its va- 
lidity to magnetic fields B for which the Kondo effect is 
destroyed, as discussed below. For sufficiently large Zee- 
man energy q/jeB, the zero-bias conductance at steady 
state is given by G = GoT{p, Vg)/2, where Go = 2e^//i 
and p is the chemical potential of the leads. The re- 
sults of Ref. [lal suggest that Eq. Q remains valid, with 
G = GoT{p, Vg), also in the absence of Zeeman splitting 
in an intermediate temperature range Ti < T < O.lw, 
where Ti is very small. Henceforth we set p = 0, corre- 
sponding to half- filled leads. 

The transmittance as a function of gate voltage is very 
small except at the poles of gno- If the ground state 
\g) has N particles for p = 0, the ground-state energy 
Eg{N) satisfies Eg{N) - eVg < Eg{N - 1) and Eg{N) + 




FIG. 2: Conductance as a function of gate voltage for the 
EAM at several magnetic fields. Circles correspond to the 
conductance of a resonant level for the majority spin only. 



eVg < Eg{N + 1). As Vg is increased, a pole in gno{Vg) 
is reached, and for larger Vg the ground state has TV -t- 1 
particles. Similarly, if Vg is lowered, a pole is reached 
when z = eVg + Eg{N - 1) - Eg{N) = 0. 

Assuming that the state of lowest energy for A^ — 1 par- 
ticles, \g{N — 1)), is not degenerate, gij may be expressed 
as a Laurent expansion around z = 0, 



9^jiz) 



+ (3ij +"ftjz + 



(4) 



where f3jj , 7jj , . . . are real coefficients and 

a, = {g{N-l)\c,MN)). (5) 

aj = |aj 1^ is the spectral weight for a local photoemission 
process. Substituting these expressions in Eqs. |(5J) and 
0, retaining terms to lowest nontrivial order in t'/r(0), 
and using that T{z) - 0(1) for z < (i')Vr(O), yields 



T{z) 



1 



4aoa„ 



o{{t'/rio)f), (6) 



1 + [w/zY (ao + a„)2 
where the half- width at half maximum peak height is 

u; =(ao + «„)(*') Vr(0). (7) 

Using G — GoT{p, Vg), this expression coincides at suffi- 
ciently low temperatures with Eq. (7) of Ref. [15j. 

If sites and n are equivalent by symmetry, then w = 
2(<')^ao/r(0) and the integrated weight / = JdzT{z) = 
TTw are both proportional to (i')^ and to ckq. If Vg is in- 
creased instead of decreased, the same result is obtained 
with Cjcr ^ c]^ and A^ - 1 ^ TV -F 1 in Eq. ©. Thus 
a single transport measurement gives simultaneous spec- 
tral information related to photoemission and to inverse 
photoemission. We stress that this information is ob- 
tained with much finer energetic resolution {peV) than 
that available by direct spectroscopic techniques (meV). 

We turn next to a calculation of the conductance for 
small or zero B, where either \g{N)) or \g{N ± 1)) is 
spin-degenerate. This degeneracy can be taken into ac- 
count accurately if w is small compared to the separation 



between groups of spin-degenerate energy levels, because 
the problem becomes equivalent to an EAM. For simplic- 
ity we take Vg such that a nondegenerate state \g{N)) is 
close in energy to the spin-degenerate state \g{N — 1)), 
but extension to other cases is straightforward. Neglect- 
ing other states, the EAM for any interacting system be- 
tween the leads is 

Ha = -tJ2 c1<tC»-i<t ~^di {t'_^C-ia + t[cia) + H.c. 
+ed ^ nda + UATid^ndi - gfiBB{nd^ -ndi), (8) 



with Uda == d^d^, t'_i = t'ao, t[ = t'a„, Sd = Eg{N) - 
Eg{N — 1) — eVg, and infinite Ua (justified because Ua = 
Eg{N) + Eg{N - 2) - 2Eg{N - 1) > w). 
The conductance is given by 

where PduU^j) is the spectral density of the effective d„ 
electrons [l7|. The conductance of the EAM computed 
as a function of gate voltage for several values of B using 
slave bosons in the mean-field approximation (SBMFA) 
[19| is shown in Fig. 2. For B — Q, G increases abruptly 
from to Go when 2; ^ 0, and remains nearly perfect for 
higher Vg due to pinning of the Kondo peak in pda{^) 
at the Fermi level. However, because the Kondo energy 
scale Tji decreases exponentially with Vg, even for very 
small B the plateau becomes a broad peak, terminated 
when Tr- < gfigB, beyond which G falls to zero. The 
abruptness of the fall is an artifact of the SBMFA, but 
the total width of the feature is well described. For larger 
B the peak width decreases, tending to w to recover the 
previous result: for gfJ-sB ^ w, G — GqT{z)/2, with 
T(z) given by Eq. ©. 

We now apply these results to the topological charge 
transition in the IHM, which is defined by 

^IHM = -ii,5](ct^i,c,;,e**/^-fH.c.)-iA^(-l)X 
-gPsB^in,^ -nii) + U^ni^nii, (10) 



where Uia- = cj^Ci^, n^ = I]^njo-, and N = Yl^^n,. We 
calculate TiVg) by exact diagonalization for the isolated 
ring with L sites and N = L electrons, with leads at- 
tached to two sites of the same energy (—A/2) in the 
presence of a magnetic flux (j) (Fig. 1), with $ = 2tt(J)/(J)q. 
The Green functions gij{Vg) are obtained numerically 
and substituted in Eqs. lO and Q. We choose L = 8, 
but similar results are obtained for any L = An {L = 
An + 2 with $ -^ <f> + tt) . The qualitative features are 
independent of the lead position n, with the exception of 
n = L/2 where the transmittance at $ = tt vanishes for 
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FIG. 3: Transmittance as a function of gate voltage for an 8- 
site IHM ring with leads at sites and 2 for tn = 1, t' = 0.2, 
A = —3, and values of U below and above Uc{L = 8) = 4.352. 



symmetry reasons. Because i?iHM is invariant under si- 
multaneous particle-hole transformation and sign change 
of A, we restrict our analysis to eVg < 0. We set t^ = 1 
as the unit of energy unless otherwise stated. 

The topological transition is present at any value of 
t/j. As a basis for our study we consider a realistic ring 
structure containing 8 QDs with only their lowest levels 
singly occupied, the centers separated by 200 nm j5| , and 
the parameters ?7, A ~ 1 meV and tfi ^ C//4 (intermedi- 
ate coupling strength). We first assume the presence of 
a strong magnetic field (i? ~1 T for a QD array), in the 
plane of the ring in order not to alter the threading flux, 
which destroys the Kondo effect. Fig. 3 shows the first 
peak in G{Vg) at A = —3 and for two values of U. For 
L = 8 and A = ±3, the topological transition occurs at 
Uc = 4.352. The peak width is parameter-dependent: for 
U = 5, w is approximately constant with decreasing flux, 
whereas for U = A, w decreases and the peak disappears 
at $ = 0. liU < Uc, the ground state \giL)) for $ = is 
even under reflection through any site (corresponding to 
a BI), while for U > Uc \g{L)) is odd (Ml or SDI) |^. For 
A = 0, the lowest-energy hole enters the system with the 
Fermi wave vector, ±7r/2, leaving an orbitally degenerate 
ground state with L — 1 particles, \g{L — 1)). For A ^ 0, 
this degeneracy is broken and \g{L — l)) is odd under re- 
flection through any even site if A is negative jig ■ As a 
consequence, for U < Uc the matrix element gq [Eq. jSJ] 
vanishes by symmetry, whence G is negligible at $ = 0. 
The flux acts as a symmetry-breaking fleld, which allows 
one to follow the flrst peak in a continuous manner until 
it disappears as $ ^ 0. 

To demonstrate that these essential features are not af- 
fected by the presence of a Kondo resonance we consider 
the 8-site ring with parameters (Fig. 4) similar to one ex- 
perimental realization [5| and B normal to the ring plane 
so that the threading flux and Zeeman splitting have the 
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FIG. 4: Transmittance as a function of Vg for a ring of 8 QDs 
described by the IHM. The parameters chosen are U — 1 
meV, A = -3C//4, t = t//4, t' /t = 0.2, and gus = 0.025 
meV/T. The curves have been shifted in Vg (cf. lower panel 
Fig. 3) such that their maxima coincide. 



same origin. Fig. 4 shows the conductance in this regime, 
calculated using the EAM in the SBMFA. The disappear- 
ance of the feature with decreasing flux remains clear. 

The topological transition may now be characterized 
using ttQ. As shown in Fig. 5 for A = 3, ao as a func- 
tion of U (cf. peak widths in Fig. 3 [ig) vanishes dis- 
continuously at t/c for $ = 0, indicating unambiguously 
the charge transition. The analogous result obtained for 
fixed U by varying A is of direct experimental interest, 
because A can be controlled by a difference in gate volt- 
age applied between even and odd sites. A finite flux 
acts as a parity-breaking perturbation and smooths the 
transition. This is the situation for artificial arrays of 
QDs, in which perfect structural symmetry is difficult to 
attain. We note that curves for all flux values cross ap- 
proximately at the same point: transport measurements 
under different applied fields can therefore help to locate 
the transition, even in the absence of perfect reflection 
symmetry. However, for a small molecule, which by defi- 
nition exhibits the $ ^ limit, the only source of asym- 
metry is the leads. 

In summary, we have considered the transport proper- 
ties of a ring-shaped interacting system connected weakly 
to conducting leads. The conductance is expressed in 
terms of the spectral density of an EAM, which illustrates 
the breaking of the Kondo effect by an applied field. Out- 
side the Kondo regime we have derived a transmittance 
formula including nontrivial interference effects. The 
conductance peaks may be characterized by their total 
weight, offering spectral information with far higher res- 
olution than conventional spectroscopic measurements. 
As an application of this result, we have demonstrated 
that the conductance can be used to detect the charge 
transition of molecules or quantum dot rings described by 
the IHM. The method is relevant in the general context 
of systems presenting phase transitions which involve a 
symmetry change of the ground state. 
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FIG. 5: Matrix element ao as a function of U for A = 3 
{Uc = 4.352) over a range of values of the flux "I". 
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